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We compute the low-temperature behavior of the specific heat of normal (non-color- 
superconducting) degenerate quark matter as well as that of an ultradegenerate electron 
gas. Long-range magnetic interactions lead to non-Fermi-liquid behavior with an anomalous 
^j ' leading TlnT^^ term. Depending on the thermodynamic potential used as starting point, 

^^ I this effect appears as a consequence of the logarithmic singularity in the fermion self-energy 

at the Fermi surface or directly as a contribution from the only weakly screened quasistatic 
magnetic gauge bosons. We show that a calculation of Boyanovsky and de Vega claiming the 
^fJ' absence of a leading TlnT^^ term missed it by omitting vector boson contributions to the 

internal energy. Using a formulation which collects all nonanalytic contributions in bosonic 
ring diagrams, we systematically calculate corrections beyond the well-known leading-log 
approximation. The higher-order terms of the low-temperature expansion turn out to also 
CO ■ involve fractional powers t(3+2")/3 qj^(\ t^q explicitly determine their coefficients up to and 

including order T'^^^ as well as the subsequent logarithmically enhanced term T'^ln(c/T). 
We derive also a hard-dense-loop resummed expression which contains the infinite series of 
r~^ I anomalous terms to leading order in the coupling and which we evaluate numerically. At low 

00 ■ temperatures, the resulting deviation of the specific heat from its value in naive perturbation 

^^ I theory is significant in the case of strongly coupled normal quark matter and thus of potential 

jir ■ relevance for the cooling rates of (proto-)neutron stars with a quark matter component. 

o 

^- I. INTRODUCTION 

Oh! 

(— I I It is well-known that long-range magnetic interactions in a degenerate electron gas lead to 

non-Fermi-liquid behavior which manifests itself in the appearance of an anomalous contribution 
- to the low-temperature limit of entropy and specific heat proportional to aT In T~^ as discovered 

r> I by Holstein, Norton, and Pincus |l| over thirty years ago. While this effect is perhaps too small 

C^ ' for experimental detection in nonrelativistic situations, it drew renewed theoretical attention more 

recently |3, y, |j| after the detection of non-Fermi-liquid behavior in the normal state of high- 
temperature superconductors [a| and in other systems of strongly correlated electrons, which may 
be due to effective gauge field dynamics (see also p, [D, |8| ) . 

In deconfined degenerate quark matter, the analogous effect can more easily be important 
because the larger coupling constant «« together with the relatively large number of gauge bosons 
increases the numerical value of the effect by orders of magnitude. In contrast to the case of a high- 
temperature quark-gluon plasma, chromomagnetostatic fields are expected to remain unscreened in 
the low-temperature limit |9| and thus lead to the same singularities in the fermion self-energy that 
are responsible for the breakdown of the Fermi-liquid description in the nonrelativistic electron gas 
considered in |l|. 

An important consequence of such non-Fermi-liquid behavior in quantum chromodynamics 
(QCD) is a reduction of the magnitude of the gap in color superconductors |9|, |lO|, |ll| which 
on the basis of weak-coupling calculations are estimated to have a critical temperature in the 
range between 6 and 60 MeV [1^. Quark matter above this temperature, and unpaired quark 
matter component also below it, has long-range chromomagnetic interactions that should lead 
to an anomalous specific heat with possible relevance for the cooling of young neutron stars as 
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pointed out by Boyanovsky and de Vega [1^, llJI- However, in Ref. [1J| these authors claimed 
that the aTlnT"^ term in the specific heat as reported in [l|, |2, y] would not exist, neither in 
QCD nor in QED. Instead they obtained a aT^lnT correction to the leading ideal-gas behavior, 
which by renormalization-group arguments was resummed into a T^+0(a} correction as the leading 
non-Fermi-liquid efi^ect on the specific heat.^ At low temperatures, such a contribution would be 
rather negligible compared to standard perturbative corrections to the ideal-gas result oc T. 

In a numerical study of the exactly solvable large-flavor-number limit of QCD and QED jl^ at 
nonzero chemical potential [la], two of us however found that the entropy at low temperature has 
a behavior suggestive of a aT In T~^ term. In Ref. [l3|, the three of us have recently reproduced 
the known aT In T"^ term in entropy and specific heat, together with further anomalous higher- 
order corrections, in an analytical calculation that should apply equally to the case of finite flavor 
number. This calculation is however organized in a form which does not allow one to compare 
directly with the calculation of Ref. [ij] where all aT In T~^ terms appeared to cancel. 

In this paper we shall therefore investigate the approach of Ref. [ij] , which derived the speciflc 
heat from a formula for the internal energy, and compare with two somewhat more direct calcula- 
tions, one using a self-consistent formula for the entropy and another using an expression for the 
thermodynamic potential that becomes exact in the limit of large flavor number. 

As we shall demonstrate, all these approaches agree eventually and do give a leading aT In T^^ 
term for the speciflc heat. In the calculation using a self-consistent formula for the entropy fSect.HHl 
the aTlnT"^ term arises as a contribution from the spectral density of the fermions with their 
logarithmic singularity in the self-energy. There are also aTlnT^^ contributions from the gauge 
boson sector, but these cancel in the end, which thus validates the (in our opinion not unquestion- 
able) starting point of the original calculation by Holstein et al. []|. On the other hand, in the 
calculation of the speciflc heat from the internal energy (Sect. IIII|) we flnd that keeping only the 
fermionic contributions leads to a cancellation of the leading qT In T^^ term, just as observed in 
Ref. \l4\. However, it turns out that in this approach the contribution of the gauge bosons to the 
speciflc heat cannot be neglected, but now contains the complete leading logarithm. 

In Sect. lIVl we describe the details of a calculation which allows us to systematically go beyond 
the leading-log approximation. Besides completing the argument of the leading logarithm, we 
find fractional powers 7"(3+2n)/3 g^^^^ ^g determine their coefficients up to and including order 
r^'^ as well as the subsequent logarithmically enhanced term T^ln(c/r). This low-temperature 
expansion requires that the temperature is much smaller than the scale set by the Debye mass. 
At temperatures of the order of the Debye mass or larger, but still much smaller than the quark 
chemical potential, a complete leading-order result which contains the infinite series of anomalous 
terms is obtained in Sect. It involves a hard-dense-loop resummed one-loop expression, which 
we evaluate numerically in Sect. IVII This allows us to study the quality of the low-temperature 
expansion, and to compare with the exact results for the large-flavor-number limit. For the sake 
of this comparison we shall throughout use the notation 



ds 



^,QCD, 



(1) 

g'^Nf, QED. 

with g the coupling constant and Nj the number of quark (or electron) flavors. At flnite Nf we 
finally evaluate our results numerically for a range of coupling which may be relevant for (normal) 
quark matter in (proto-) neutron stars, with the finding that there is an interesting range of 
temperature where the anomalous specific heat exceeds significantly the ideal-gas value. 



Resummation of the aTlnT ^ term along the hnes of Ref. |l4J| would have led to a yi+Oia) term instead. 



II. ANOMALOUS SPECIFIC HEAT FROM THE ENTROPY 

The specific heat C^ per unit volume is defined as the logarithmic derivative of the entropy 
density with respect to temperature at constant volume and number density 

This is related to derivatives of the thermodynamic potential with respect to T and /i by [1 



but at low temperatures one has 

C. = t(|^) +0{T% (4) 

so that both Cv and S contain the same TlnT term, if any. 

The entropy as first derivative of the thermodynamic potential with respect to T is in some 
important respects a simpler quantity than Q. In gauge theories with fermions the latter is g iven 
by the following functional of the full propagators D (for gauge bosons), S (for fermions) J19t |. 

(3n[D,S] = ^TT\nD-^ - ^Ti HD -Ti In S'^ + TrJ:S + ^[D,S], (5) 

where <i> is a series of 2-particle-irreducible (skeleton) diagrams and where for simplicity we assumed 
a ghost-free gauge. 

Using the fact that r2[D, 5] is stationary with respect to variations of D and S, one can derive 
an expression for the entropy which to two-loop order in the skeleton expansion is entirely given by 
propagators and self-energies 20|,|2l|. Neglecting the longitudinal gluon mode, and the antiparticle 



contributions in the fermionic sector, this reads 



^ /apx f ,'K 



\9TJ^ J (2») 



INg }^ ' (Im InZ)^^ -ImnTReL>T) 



+ANNf ^ll^' (Im In 57^ + Im S+ Re S^ 
oT 



+ S\ (6) 



where D^^ = -w^ + fc^ -^Ut, 5+^ = -uj + k-T.+, nb{uj) = (e'^/^ - l)-\ and n/(a;) = {e^'^-t^)/'^ + 
l)-i. 

In the original derivation of the anomalous specific heat in QED by Holstein et al. yj, only the 
term involving Im In 5"^ in the quark part had been taken into account, by way of reference to 
a formula by Luttinger |22| (Eq. (46) therein), which is fully justified actually only for standard 
Fermi-liquid systems. 

Since the calculation of Ref. jlj has been questioned by Boyanovsky and de Vega Jl4^, we 
consider the more general expression above. For the moment we are interested only in terms of 
order ^f^TlnT, so we can neglect 5', which vanishes at two-loop order in the skeleton expansion 
[2i| and should therefore only give contributions which are suppressed by an additional factor of 



g'. 



A. Quark part 

In Eq. Q we have the fohowing contribution from the quarks, 

d^K dnf[uj) 



'(q) 



ANN 



^ J (27r)4 dT 
Im hi {—uj + k — E-I-) + IniS+ Re ■ 



-ijj + k — Ti_ 



■K 



^-^NNf I dkk' 



duo 



dnf{uj) 
dT 



Im In (—u) + k) — Re S+ Im ■ 



1 



-co + k 



(7) 



where we have performed an expansion with respect to E+, keeping only the free term and the term 
corresponding to a single quark self energy insertion^ . The free term gives the particle contribution 
to the free fermionic entropy density, 



5S^-iViV/^. 



(8) 



In the last term in Eq. ((T)) the factor Im 1/(—uj + k) forces the self energy to be on the mass shell. 
Using the expression for S_|_ given in 10|, |23| , 



S+=^ 



9'Cf 



{lo — fi) In 



M' 






247r2' "^^ \{uj- n)^J ' ' 127r 
which is nonanalytic in uj (but not with respect to k J_lfl|), we obtain 



(^) "^io dT 247r2 



{k — fi) In 



M' 



(k-fi) 



2' 



(9) 



(10) 



With the substitution k = Tz + /i we find that the integral is dominated by small values of z, and 
therefore we may send the lower integration limit to — oo. Then we obtain at order TlnT 






(11) 



This result agrees with the one of Holstein et al. [jj| after correcting a factor of 4 therein, 
done previously in Ref. 



as 



B. Gluon part 



The gluon part 5(g) is given by the first line of Eq. 0. Using the relation 

/ Inill \ 
Im InD-^ = arctan 7--— r - 7re(a;)(9 (-ReL»"^) , 
\ Re D ^ / 



(12) 



^ Note that, diagrammatically, the part with a single self energy insertion corresponds to the gluon ring diagram of 
section liiBI 



we write 5(g) = S^cnt) + ^(n) + 5(poie) , with 

5(cut) = 2iV, y J^.^j- arctan (^ ^,_^,_R^nJ ' ^^^^ 

-(n)^--./|^^-n.Ke-,,3A__, (M) 

(^^T^^e(^)^ l'^' - fc' - Renr) , (15) 

where we again neglect the contribution of the Debye-screened longitudinal gluons. For the cut 
term we use the approximation a; <C A:, because it can be checked that including terms of higher 
order in u) would only produce terms of higher order than Tin T (see section IIV|) . In this region 
we have 

Introducing an UV-cutoff kc for the moment, we obtain 

^^-'^ "2^/0 ^ loo ^"^T- "^^'"" [-^^) ■ (^^) 

In order to evaluate this integral we make the substitution y = to/T, x = Airk^ / {g^^^'^T) . Keeping 
only the term of order TlnT, we obtain the cutoff-independent result 

The determination of the constant M' requires a more accurate calculation and will be carried out 
in Sect. |TV| 

Next we evaluate 5(n) ■ Following similar steps as in the computation of 5(cut) i we find 

We observe that at order TlnT this expression just cancels the contribution from Eq. (|18() . 

Finally we consider the pole part. To leading order, i.e. in the hard dense loop (HDL) approx- 
imation 24l. I25I. l2q. I27 1 ■ we have at low temperature fidHx/dfi ~ 2112^, and therefore 

= -i^N, I ^5^(^2 _ j2)^,^) J (^2 _ j2 _ K^n^) , (20) 

where we have discarded contributions ~ T^ which are negligible in the low-temperature limit. 
Using 13 

e{uj) 5 (Re T»^^) = Zrik) [ 5{uj - uoxik)) - 5{u + urik)) ] (21) 



we find 



.^-^(poie) _ 2iV. r...2dn,{u:Tik)),, ,,,, ,2, 



We can estimate tliis integral as follows. Assuming T <^ ujp oc gcSfJ-, we have the inequalities 
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JO 



and 
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2 Jo dT 2^^ dT 8r2 ^ ^ 
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dke °^^^^ " u:T{kYZT{k)< j dkk'^^^^^iiOp + ky^^ 



<U^'dkk{u;, + kf^^^ + U^dkk{co, + kf 



9T 



17a;, 
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24r2 



£„-a;p/T 



e-"'^/' . (24) 



Apart from terms ~ T^ which are dropped in the derivative with respect to /i, this crude esti- 
mate (which we shall refine in Sect. ^ below) shows that the pole contribution is exponentially 
suppressed, essentially because of ujt > ^p- 

C. Result 

In total we find the following result for the entropy at low temperature, 

5..„.5,„.A.iV,^.?i^I„(^). (25) 

From Eq. ()18|) we see that the TlnT term can also be obtained by starting only from the 
expression 

S ^ NNf^ - 2N, I -^^Im \nD;,\ (26) 

^3 V (27r)4 dT ^ ' ^ ' 

with D the resummed gluon propagator. This formula corresponds to integrating out the fermions, 
as has indeed been done in the approach of Ref . [4] . 

On the other hand, we see from Eqs. ((TJ and (jlip that one equally gets the correct result by 
using only the purely fermionic expression 

(27r)4 "/t ^ "^ ^''^"'^ + ^ " ^^^+^' ^^^^ 

which justifies the starting point of Refs. [l|. 



III. SPECIFIC HEAT FROM THE ENERGY DENSITY 

In Ref. [ij], which did not find a term g^TlnT^^ in the specific heat, the starting point was 
instead the internal energy density. 

The energy density can be obtained from the expectation value of the energy momentum tensor, 

U = ^ld^x{T''{x)), (28) 

and the specific heat is then given by 

^-(^). '-) 

Here the temperature derivative has to be taken at constant particle number density, in contrast 
with the calculation of the low temperature specific heat in the previous section, where all tem- 
perature derivatives were taken at constant chemical potential, see Eq. 0. In [ij] this fact was 
mentioned as explanation for the disagreement with the previous calculation, but, as Eq. Q) makes 
clear, this could only affect terms of order T^ in the low-temperature expansion. Indeed, we shall 
show now that a complete calculation based on the internal energy also leads to a g'^TlnT~^ in 
the specific heat. 

In [ij] the specific heat is computed using the following formula for the total energy density. 






U = 2 I duj I -i——^nf{uj)ujp+{uj,k), (30) 



where p+ is the spectral density of the positive energy component of the quark propagator (see 
below). It should be noted that this formula is incorrect even for a theory with only instantaneous 
interactions of the type 

Hint = 2 / (i^xd^x'lpl{-X.t)lpj^{xt)Vaa',l3l3'{'X - x')V'/3' (x'i)'0a' (^t), (31) 

in which case the correct formula reads |2 

/r ^3^ 2 

du^ J J^M^) 2^^ + ^)-o+(^' ^)- (32) 

The anomalous behaviour of the specific heat comes from dynamically screened interactions, whose 
non-instantaneous character cannot be neglected. It might be rather difficult to generalize Eq. H32|) 
directly for non-instantaneous interactions, because one would have to use an effective Hamiltonian 
which is nonlocal in time. Therefore, we will use the full energy momentum tensor of QCD without 
integrating out the gluons. 

The energy momentum tensor can be written as a sum of three distinct pieces, 

(q) (g) (int.)' ^ '' 

corresponding to the quark part, the gluon part, and the interaction part. The contributions of 
these parts will be evaluated in the following subsections. We will neglect gluon self interactions 
and ghost contributions, since they give only higher order corrections at low temperatures. 



A. Quark part 

The quark part is given by 

/ 

where we have written exphcitly the sum over flav or sp ace. This is the (only) contribution which 
is taken into account by Boyanovsky and de Vega [lj|. We now repeat their calculation, but for 
simplicity without the renormalization group improvement of the quark propagator proposed in 
|l4| . Taking into account only the positive energy component of the quark propagator, we flnd 

dioj —-^nf{Lo)ijp+iu;,k), (35) 

where the spectral density is defined as p+ = ilmS'+. 

In order to obtain the specific heat from Eq. H29() . we have to determine first the temperature 
dependence of the chemical potential from the condition 

dj\f _dM dfidM 



^f = 2NNf I duj I j^^^nfioj) p+ico, k) (37) 



where the particle number density J\f is given by 

d^k 

{2-k)" 

(up to anti-particle contributions.) We expand N with respect to g, 

M = Mo + glsM2 + ... (38) 

The free contribution A/q is given by 

In A/'2 we are only interested in contributions which contain ln(M/T). Such terms arise from 
infrared singularities caused by the transverse gluon propagator, which are dynamically screened. 
This corresponds to scattering processes of quarks which are close to the Fermi surface. Therefore 
the anomalous terms come from the region A; ~ a; ~ //, where S^- is given by @. Subtracting the 
temperature independent part, we find then 



gl^M2 = 2NNf I dw I -^ (n/(a;) - e{p - uj)) 6{uj - k + Re S+) 



(40) 



where we have approximated the spectral density by a delta function, since the imaginary part of 
S+ turns out to be negligible compared to its real part. The integration can be performed easily, 
with the result 

We notice that this result is consistent with the result for the entropy, Eq. (|25() . Now we can solve 
Eq. p6|l at low temperature, 



dT 3/i 18// \T J 



The approximate solution to this differential equation is given by 

Eqs. (|42|) and (|43|) correctly reproduce the beginning of the perturbative expansions of the corre- 
sponding formulae in Q (Eqs. (2.37), (2.38)). 

For the specific heat we obtain from Eq. (|35() . following the same steps as in the calculation of 
(W/dT, 

C.,„.A.A.,.^T.A.A.,^|.%!^.„(^)^ (44, 

Using Eqs. ()42() and 1)43^ we find that the TlnT-terms cancel, 

a(q)=^iViV/^+0(T3), (45) 

as stated in [lj|. We should emphasize that this cancellation has nothing to do with the nonper- 
turbative renormalization group method which is employed in [ij] (and which has recently been 
criticized in Ref. |30|). 

The authors of Ref. 1J| also determined a contribution which prior to the renormalization group 
improvement corresponds to a term of order g^T^ ln(M/T) . This type of nonanalytic terms however 
appears already in regular Fermi-liquids l^lj and moreover is subleading to ordinary perturbative 
corrections g^^^T at low T, which have not been evaluated in |14| . 

B. Gluon part 

We now turn to the gluon part of the energy density, which has explicitly been neglected in 
Ref. _14]. This is given by 

Tg = 1 (E" • E" + B'^ • B") . (46) 

Neglecting gluon self interactions, and keeping only the transverse part of the gluon propagator, 
we obtain 

^fe) =^ 2iV^ / 0^M^) I™ ((^' + ^') Dt) ■ (47) 

The pole contribution to this integral is again exponentially suppressed, therefore we only have to 
consider the cut contribution. At low temperature the temperature dependence of the gluon self 
energy can be neglected and we find 



Sp^ffM^'-' 



(2^)^7o ^ (^2_p)2+(^^) 

where we have dropped less infrared-sensitive contributions not involving the Bose distribution. 
After the substitution co = Ty, k^ = g'^gfi'^Tx/{4TT) the integral can be readily done, with the result 

^i.)^'^^§^HM/n (49) 
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which gives the following contribution to the specific heat at order TlnT, 

Ch.)^'^^^HM/T). (50) 

Again the determination of the constant M would require a more accurate calculation, similar to 
the one in section HVl 

C. Interaction part 

The interaction part is given by 

T^S.t)=9Y.^l'^lT^i^- (51) 

/ 

The expectation value of this term is essentially given by the resummed gluon ring diagram as in 
[15|. However, here only the longitudinal component of the gluon propagator appears in the loop. 
This mode is subject to Debye screening, so it can contribute only to the normal Fermi liquid part 
of the specific heat. 

D. Result 

We have thus found that the only contribution to the specific heat at order TlnT when cal- 
culated along the lines of Ref. [ij] comes from the gluon part, Eq. (|5U() . While this confirms the 
observation of Ref. [ij] that the quark contribution of order g^T In T cancels against a similar term 
in the temperature dependence of the chemical potential at fixed number density, it shows that 
the (explicit) neglect of the gluon contribution to the internal energy in [lj| is not justified. 

Curiously enough, the gluon contributions can be neglected when calculating the anomalous 
specific heat from the entropy functional ®. But the calculation of section^can also be viewed as 
receiving a net contribution only from a purely gluonic term H26|) , since the anomalous contribution 
contained in the S'(n) part is equal but opposite in sign to the fermionic part S'(q) . In this respect 
the calculation based on the internal energy is perfectly in line with the calculation based on the 
entropy. The different possibilities for organizing the calculation of the anomalous contribution to 
the specific heat thus correspond to 'integrating out' first the fermionic degrees of freedom or first 
the gluonic ones. 

IV. HIGHER ORDERS IN THE LOW-TEMPERATURE SPECIFIC HEAT 

In this section we shall evaluate higher terms in the low-temperature expansion of the specific 
heat which go beyond the leading log approximation. A convenient starting point is the follow- 
ing expression for the pressure, which can be viewed as the result of having integrated out first 
the fermionic degrees of freedom, and which thus concentrates on the effects of the only dynam- 
ically screened transverse gauge bosons in the low-temperature expansion of the thermodynamic 
potential: 

P = NNf ( :r^ + ^ \- 



127r2 6 180 

T ^2 



'''' I (27r)3 io vr 



Tih + -]Im InT)^^ - -Im InTT^ 
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FIG. 1: Integrand for the g-integration g^Im \n{q'^ 



ll 



Ht + Hvac) for A* = MMs/2 = 1, go = 0.2, g^ 



2 
eff 



1. 

The sohd hne shows the exact result that follows from the full one-loop self energy expressions at T = 0, 
the dashed line shows the result with the approximations of Eqs. ifS^ and (|S7|) . The parameter go = 0.2 is 
chosen this large as to clearly show the three different ranges. As discussed in the text, the main contribution 
only comes from region II. 



[nb + -]Imln-2-^^2 



q^ 



% 



1 D-^ 

ilm In ^'^^ 

2 g2_^2 



+ +0(<7^n+0(5V^). 



(52) 



Here the inverse gauge boson propagators are given by Dj, = q^ — q^ + Hy + Hvao Dj^ 



+ 11^ + Ilvac, and D^^ 



q^ + Ilvac) where TiT,L,vac are the matter and vacuum 



contributions to the gauge boson self-energy produced by an undressed one-loop fermion diagram. 
This expression becomes exact in the limit of large flavor number Nf and has been studied in 
Refs. |l5l . Ilq and used to test the behavior of perturbation theory at finite temperature and 
chemical potential. At finite Nj, Eq. ()52|) with 11 including also the leading contributions from 
gluon loops still collects all infrared-sensitive contributions up to and including three-loop order 
|32| . We shall however find that all contributions from gluon loops to 11 enter the specific heat 
only at and beyond order g^^T^ , and thus will be negligible for T <^ jjl compared to the terms we 
shall keep. 

In the following we will always drop temperature independent terms in the pressure, since they 
do not contribute to the specific heat at low temperature. 



Transverse contribution 



The rife-part of the contribution of the transverse gluons to the pressure is given by 



Nr, 



Sc 



(27r3) 



dqo 



IT 



rifelm ln{q 



+ Ht + n^ 



(53) 



As long as r <C /U, it is sufficient to take the self energy at zero temperature, which is given by 
a single fermion loop. From the explicit form of the self energy J33l | we see that the g-integration 
naturally splits into three regions: q < qo (1), qo < q < 2fi — qo (II), and 2iJ, — qo < q < 2/1 + qQ 
(III), see Fig. n 

Usually region I contains the ideal-gas pressure of the gluons, 7r^T^/45 per gluon, and pertur- 
bative corrections oc g'^^fjpT'^. However, at low temperatures T <C Wp ~ geflfA*) these contributions 
are suppressed by a factor e~^vl'^ which goes to zero with all derivatives vanishing and thus do not 
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contribute to the low-temperature series. All other contributions from region I are suppressed by 
further powers of ges ■ 

In region III we may expand the self energy about q = 2fi and qq = 0. Then one finds that 
this contribution is of higher order in g^s ■ We conclude that we may restrict our attention to the 
region II. 

In region II the Bose-Einstein factor and the leading term in the gluon self energy set the 
characteristic scales. Since qo ^^ T '^ geS^^, this is dominated by the well-known [2J| dynamical 
screening pole at imaginary q with \q\ ~ ((^^g/i^go/l^vr))-^'^ and we have 



qo 



^~(5eVT)V3 



(54) 



in the infrared-sensitive part of region II. 

We shall therefore perform an expansion with respect to a parameter b defined by 



b:-- 



T 
geSfJ' 



1/3 



(55) 



It turns out that the following approximation of the gluon self energy is sufficient through order 
T'^ In T in the entropy (see Fig. ^ , 



vr"^ \ q"^ 



3g4 



Att 



(56) 

(57) 



The first two terms in both lines are the leading terms of an expansion of the HDL self energy 
2J, |2a, |2a, 1231 , in powers of q^. Naively counting powers of b in the integrand one would conclude 
that only these terms are responsible for the terms of order b^ through 6^^ in the pressure (perhaps 
with additional factors of In 6). In principle this is correct, but one should keep in mind that the 
integration limits of the g- integration depend on /x and goi which might invalidate a naive power 
counting. However, it turns out that there is only one instance where a term which is formally 
suppressed in the naive power counting scheme has to be included in the self energy (see 'Sfr] for a 
rigorous proof): It is the term of order b^ in the pressure, where one also finds a contribution from 
the last term in Eq. (|57|) . which is beyond the HDL approximation (but it is still included in the 
large-A'^y limit.) The b^ term plays a special role anyway, as this is the only term where we also 
get a contribution from the non-n^-part (see below). 

Introducing dimensionless integration variables x and y via q^ = b^gcsfJ-y and q = 
bgcs fJ'ix / (Att))^'^ , we find after expanding the integrand with respect to b, 
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(58) 



with Xmin = Airb^y^ and Xmax = 47r(2 — b"^ g^gy)^ / (bgeff)^ . In the coefficients of this expansion we 
have written down only those terms which do not ultimately lead to terms that are suppressed by 
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explicit positive powers of g^fi- The integrations are now straightforward, and we find 

b^ (, ( 327r 
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(59) 



The evaluation of the constant ct is a bit more involved because one has to sum up an infinite 
series of contributions from the infrared region. This calculation is performed in the Appendix and 
leads to an integral representation of ct that we have been able to evaluate only numerically, with 
the result ct = -0.001786743 .... 

B. Longitudinal contribution 

The nf,-part of the contribution of the longitudinal gluons to the pressure is given by 
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(60) 



As in the previous section the g-integration decomposes into three parts. 

Again, contributions from region III are suppressed by explicit powers of g^s compared to those 
of the other regions. 

The dominant contribution comes from region II, q^ < q < 2^ — q^. Now the characteristic 
scales are 



go ~^, 



ffcff^, 



(61) 



because of Debye screening of the longitudinal gluons with leading-order mass niD = gcS^^/T^- In a 
similar way as in the previous section, the gluon self energy can be approximated as 
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(63) 



We introduce dimensionless integration variables y and z via q^ = b^gcsfiy, q = gcSfJ-^/T^- Then 
we find after expanding the integrand with respect to b, 
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with Zmin = fe'^yTT and 



(2/5cflf — b^y)TT. In the coefficients of this expansion we have written 



down only those terms which do not lead ultimately to terms that are suppressed by explicit 
positive powers of g^s. The integrations are now straightforward, and we find 
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(65) 
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The constant c^ can be determined by summing up IR enhanced contributions in analogy to the 
constant ct of Sect. lIV'X| and its integral representation is given in Eq. ()A8|) . Numerically, we get 
c^i Ri 0.11902569216.... 

In contrast to the case of transverse polarizations, there is now however also a contribution from 
region I, of the order of g^gfi^b^'^ = T^. The term involving Im ln(g^ — Qq + ^L + Hvac) is again 
exponentially suppressed for T <^ LOp ^ 5eff A*- However, at these temperatures the term involving 
Im ln(q^ — Qq) contributes the equivalent of an ideal-gas pressure of one bosonic degree of freedom, 
but with negative sign, leading to 



vr 
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Ng 90 
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(66) 



so that CL = el - 7r^r4/90 = 0.009363421 .... 

C. Non-rtf, contribution 

The nonanalytic terms in the low-temperature expansion of Eq. ()52() all come from the parts 
of the integrals involving the Bose distribution n;,. The non-n^ parts in Eq. ()52() are less IR 
singular and can be calculated by expanding out the self energy diagrams. We can determine their 
contribution by the observation that at two-loop order also the ni, part is IR safe and given by 
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(67) 



In this integral we have two contributions, one from the real parts and one from the imaginary 
parts of the gluon self energies. One finds that these two contributions cancel precisely at the 
order g'^giJ?T'^. (As above one finds that the T = gluon self energies are sufficient at this order.) 
Therefore the two-loop non-nj, contribution has to be equal to the standard perturbative result at 
order gl^^i^T'^ EF 
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(68) 



D. Combined result 

Our final expression for the leading temperature-dependent contribution to the interaction pres- 
sure in the regime T <C ^/x is contained in 



AP = P - PO = P|I + Pi + Pj} + Pnon-n, 

where P^ is the ideal-gas pressure, and explicitly reads 
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where the constant c is given by 



7B-90 



C'(4) 



vr^ 



31 

12 ^ 2048 



2567r2 - 367r4 + Stt" 
X <! 37r^(12 - vr^) Invr + 128(71^ - 8) ln(47r) + 3tt^{29tt^ + 32) 

1024 + vr^ - 647r2(2 + Sz^) + 127r''(-l - z"^ 



-nir-" 



dz 



Zz-^ 



247r3; 



+2z arctan 



4.099348512039 .... 



7r(l-z2) 



2z+(z2_l)ln(^f±i 



+ z arctan 



2z 1 , 
h -In 

vr vr 



z + 1 



(71) 



The terms involving logarithms and fractional powers of T all come from the cut contribution of 
region II, whereas P^ — P^ = —^P^ + 0{e~''^''^) only contributes to c. 

From Eg. (|7Uj) one can obtain the entropy density through S = {dP/dT) and the specific heat 
through 
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with N = {dP/d^)rp, which in the ideal-gas limit reads 
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(73) 



For the interaction part of the specific heat only the logarithmic derivative of the entropy in 
formula ((7^ contributes, and is given explicitly by 
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(74) 



We remark that in Eq. (|7i|) one could replace /i with /i(T), as given in Eq. (|l3|) . as this would 
modify the result only beyond the terms of order T"^ since T/{gfi) <C 1. 

Using the method described in Sect . IIV'XI - IIV CI one can in principle compute the coefficients 
of higher terms in the expansion of C^ with respect to b. This is straightforward for the coefficients 
of the fractional powers and the logarithmic terms, where one only has to include higher orders 
in the expansion of the HDL gluon self energy with respect to qo, see Eqs. (|^l)|) . (|^7|) . (|H^ . (|HH|) . 
For the terms of order T^, T'^, T^ etc., however, one has to sum up again IR enhanced terms, in a 
similar way as in the calculation of c described in the Appendix. 

The low-temperature expansion that we have carried out has assumed that T <C 5eflF^ as well 
as QeS <C 1. If we set T/ij, ~ g^^ with 5 > 0, we find that the terms in the expansion (|7Hl 



correspond to the orders g^^ ln{c/ges), ffeff 



3+(5/3)(5 3+(7/3)<5 



g^s , and g^^ ln(c/5reff), respectively, with 
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a truncation error of the order g^^ from higher terms in the ring diagrams. We have neglected 

perturbative corrections to these terms, which at a minimum arise at the order g^Q ln{c/ges)- 

One might suspect that higher order terms could involve also higher powers of g'^glnT, which 
could resum into a leading term lU^T^+'^'-^cftO. However, it has been argued in |2l] that the leading 
^fg^rinT"^ is not modified by higher order corrections in QED, and this has been corroborated 
recently by the authors of Ref. |30|] using a high-density effective field theory also applicable to 
QCD. It can therefore be expected that the leading term in the above low-temperature series 
remains valid even when the temperature is so low that g^f^ ln{gesfi/T) S> 1. 

On the other hand, the higher terms of the low-temperature expansion involving fractional pow- 
ers 7"(2n+3)/3 -^[^i^ n > 1 remain more important than the undetermined perturbative corrections 
(which are suppressed by explicit powers of g'^g) only when 6 < 3/n. 

V. HDL RESUMMATION 

As we have seen in Sect. IIVI the nonanalytic terms in the low-temperature expansion of the 
thermodynamic potential are determined by HDL contributions to the gluon self-energy. Terms 
beyond the HDL approximation are relevant for contributions from hard momenta g ~ /i, yielding 
a term of order g^^fi'^T'^ in the temperature-dependent part of the pressure. However, this is a 
perturbative piece that can be identified as a two-loop contribution without the need for resum- 
mations. When this contribution is subtracted from the full one-loop expression, the remainder is 
dominated by soft momenta q <^ fi and the HDL approximation is sufficient. 

In this section we shall consider the full HDL-resummed ring diagrams, which allows us to relax 
the requirement T <C geSfJ-, under which the above low-temperature series is meaningful, to only 
T <^ fi. When expanded around T = 0, the one-loop HDL-resummed thermodynamic potential 
contains all the anomalous terms of the low-temperature series (|7()|) . However, as we have already 
seen there are also terms from region I which behave as ~ e~'^p''^ and thus do not show up at 
any finite order of the low-temperature series. Nevertheless, such terms become important for 
T ~ gcflf^- By numerically evaluating the full HDL-resummed expression we can capture their 
effect as well and study the behavior of entropy and specific heat for all temperatures T <C /i. 

A. Separation of hard and soft contributions 

In the transverse sector, the one contribution in Eq. (|58|) from a non-HDL term in the gluon 
self-energy can also be written as 
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Similarly, in the longitudinal sector the non-HDL contribution to Eq. H64() is 
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and again Inill]^ — > for q ^ qQ. 

The HDL part of the gluon self-energy, explicitly given by 
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(80) 



with mo = QeSfJ'/'^ is the leading-order contribution at small qo,q <^T provided T <^ fi. In order 
to retain all contributions that are nonanalytic in T at T = 0, the HDL self-energies need to be 
kept unexpanded in 
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where we have dropped contributions from region HI as being suppressed by explicit powers of (7cff- 
Individually, the transverse and the longitudinal contributions depend logarithmically on the 
upper integration boundary, because 
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at large q. However, the combined expression (|5T|) is saturated by soft momenta g <C /i, and the 
upper integration limit 2// can be sent to infinity. This just amounts to dropping terms that are 
suppressed by explicit powers of ^eflf- 

The only other contribution that needs to be taken into account is Pnon-nj,, which as discussed 
in Sect. IIVCI can be treated perturbatively to the order under consideration. Put together, the 
final result for AP = P - P° is 
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(84) 



where we have subtracted the ideal-gas pressure of the gluons by including a free propagator in 
the argument of the first logarithm. 
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As indicated, this expression provides the leading terms in the temperature dependent part 
of the pressure and therefore the leading terms in entropy and specific heat. For T <C gcSfJ- the 
contribution to the pressure (as opposed to entropy and specific heat) is subleading compared 
to the three-loop result for the zero-temperature pressure obtained by Preedman and McLerran 
|37|, Baluni '3e], and Vuorinen [3j|. However, when T > gcSfJ-, its magnitude is comparable to 
g^ff^^ln{c/g), the highest known term in the perturbative result of the T = pressure, and thus 
provides an extension of the result of Refs. |32l . l37l . l38| to non-zero temperatures in the domain 



B. HDL quasiparticle pole contribution 

In region I, i.e. < q < qo, the HDL propagator has single poles for qo > ojp = mj;)/\/^ at real 
qo = uJT,L{q), which allows us to carry out the qo integration, yielding 
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At small temperatures T <C Wp, the contribution from transverse polarizations is suppressed by 
a factor e~^pi'^ and thus does not show up at any finite order of a low-temperature expansion in 
terms of powers and logarithms of T. Using that for small q-^ojp the dispersion law of transverse 
gluons is given by |3J] 



one can calculate the leading term as 
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As discussed already in Sect. lIV'Bl the low-temperature contribution of the longitudinal gluons 
involves a contribution oc T^ as well as exponentially suppressed terms. Using that 
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At temperatures ujp <^ T <^ fi the dispersion relation of the longitudinal gluons approaches 
the light-cone exponentially, which gives an equally exponentially vanishing contribution to the 
pressure. The transverse gluons, on the other hand, tend to the mass hyperboloid wt(<?) ^• 
y^g^ _|_ jj^2^ with rn^ = rnj^/l, yielding 
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C. HDL cut contribution 

For q > qo, the HDL self-energies have an imaginary part oc rrijj = g^^ij?' / ir"^ which corresponds 
to cohisionless Landau damping of hard fermions. At low go ~ ^ ^ 5eff/^, this provides the 
dynamical screening of quasistatic magnetic fields which is responsible for the anomalous terms in 
the low-temperature expansion of the specific heat, Eq. H54|) . 

At higher temperatures Wp <C T <C /i, it is instead electric Debye screening which gives the 
dominant contribution from soft momentum scales 
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(91) 



(At still higher temperatures T ^ fi, one eventually has to replace the value of the HDL Debye 
mass by the hard-thermal- loop result 28]). 

The appearance of a contribution oc m^ is traditionally referred to as plasmon effect. Our full 
HDL-resummed result ()84() gives a unified description of this longitudinal plasmon effect with the 
anomalous (non- Fermi-liquid) contributions from transverse quasistatic fields which are only weakly 
screened, and interpolates between these two different effects. As the temperature is lowered, the 
longitudinal plasmon term m^T which is linear in T gradually disappears and gets replaced by a 
quadratic term oc m'j-)T'^lii{mD/fJ-), cf. Eq. ()65() . This then combines with the leading anomalous 



term oc m'j^T'^ ln(T 



1/3 2/3 
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/fi) from the transverse sector, cf. Eq. (|59|) where rriD enters through 



dynamical screening, and whose logarithmic dependence on the hard scale fi cancels that of the 
longitudinal sector. 

VI. NUMERICAL RESULTS 
A. Full HDL result versus low-temperature expansion 



We shall now turn to a numerical evaluation of entropy and specific heat following from the HDL- 
resummed pressure (|84j) . The corresponding expression for the entropy density is given explicitly 
by 
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(92) 



where S^ is the ideal-gas entropy density. Eq. ()92() represents the leading interaction term at weak 
coupling for all T <^ fi. It is essentially given by one universal function of the dimensionless variable 
T/{gcsiJ-), which we define through 
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(93) 



and which we have normalized such that the ordinary perturbative two-loop result 3o\ for the 
low-temperature entropy corresponds to & = —1. 

In Fig.|21we display the individual contributions to © as provided by transverse and longitudinal 
quasiparticle poles (region I), and the Landau damping cut (region II). As one can see, the latter is 
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FIG. 2: The function 6{T/{geSiij) which determines the leading-order interaction contribution to the low- 
temperature entropy. The normalization is such that S = — 1 corresponds to the result of ordinary perturba- 
tion theory. The dash-dotted line shows the contribution from region II, comprising HDL Landau damping 
and hard contributions; the two dashed lines give the transverse (T) and longitudinal (L) quasiparticle pole 
contributions of region I. 

responsible for the anomalous behavior of an interaction contribution which is positive for T/{gfj,) < 
0.0404 . . . and is well reproduced by the low-temperature series involving logarithms and fractional 
powers of T/{gggii). The quasiparticle contributions, on the other hand, behave as T^ at small 
temperatures, but cannot be further expanded about T = due to terms involving e~'^^''^ . 

We have in fact been able to perform the required numerical integrations with sufficient accuracy 
to explicitly check all the coefficients of the low-temperature expansion calculated in Sect. lIVl (and 
further ones up to order T^ In T [3a| ) • 

In Fig. 131 we compare the first few orders of the low-temperature series with the full HDL result. 
The low-temperature result up to and including the T^ In T contribution to the entropy is a good 
approximation for T/{gegfi) ^ 0.04 where the anomalous contributions dominate; for larger T the 
non-expandable e~'^''''^ terms in the quasiparticle pole contributions become important numerically. 

B. Comparison with nonperturbative large-Nf results 

When applying our results to QED as well as QCD, the range of T//_i where one finds an excess 
of entropy and specific heat over the ideal-gas result will be the larger the higher the coupling is. 
However, we then have to expect more important perturbative corrections which are suppressed 
parametrically by further powers of g^. In order to assess their importance, we compare with the 
special but exactly solvable case of infinite flavor number. This is done in Fig. ^ for the three 
valuesfleff (/"MS = 2/u) = 1,2,3 where the heavy dots give the nonperturbative large-A^j result of 
Ref. [lal and the full line represents the full HDL result (solid line, denoted by HDL"*" to remind 
of the inclusion of hard, perturbative terms). Also given is the low-temperature series up to and 
including the T^ In T contributions. 

While this certainly does not allow one to predict the reliability of our HDL result for real, 
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FIG. 3: The first few orders of the low-temperature series for the entropy as determined by Eq. (|70|l in 
comparison with the full HDL-resummed result. 



0.003 



0.002- 



0.001- 



-0.001- 



-0.002- 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 



T 



FIG. 4: Complete entropy density in the large- iV/ limit for the three values gcsiP-MS = 2^) = 1,2,3 (heavy 
dots), compared with the full HDL result (solid line) and the low-temperature series up to and including 
the T'^lnT contributions. 



fimte-Nf QCD, it should give an idea of the errors to expect at least. Interestingly enough, in 
large-A'^y QCD the higher-order corrections seem to increase somewhat the range in T//u where 
there is an excess of the entropy over its ideal-gas value. 

By the same token, in Fig. [S] we display the renormalization scale dependence of the HDL result 
by varying the renormalization point by a factor of 2 around a central value of /^ms = 2//. 
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FIG. 5: Complete entropy density in the large- A''/ limit for the three values gcsifius = 2^) = 1,2,3 (heavy 
dots), compared with the HDL-resummed result when in the latter the rcnormalization scale is varied by a 
factor of 2 around /Ums — 2/i. 







QCD strict 
QCD th.dyn.cons. 
QED strict (x20 enlarged) 
QED th.dyn.cons. (x20) 
Pert.th. strict 



0.12 0.14 



T 

/i 



FIG. 6: The HDL-resummed result for the specific heat Cy, normalized to the ideal-gas value for g^s — 2 and 
3 corresponding to as ~ 0.32 and 0.72 in two-flavor QCD, and gcff ~ 0.303 for QED. The results labelled 
"strict" does not include anomalous contributions in the second term of Eq. H72|l where the would be of 
higher order in geff, whereas "th.dyn.cons." refers to a less systematic but thermodynamically consistent 
evaluation. The deviation of the QED result from the ideal-gas value is enlarged by a factor of 20 to make 
it more visible. 



C. Specific heat 

In Fig. we finally evaluate our result for the low-temperature specific heat at constant baryon 
density, C^. As can be seen from Eq. (|72() . this is a nonlinear functional of the thermodynamic 
potential. However, for /i ^ T and to leading order in (/eff the anomalous contributions provided 
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by Eq. ()84() enter only through the logarithmic derivative of the entropy, and the nonlinear terms 
in Eq. (|72|) need only include the ideal-gas result. This defines the results in Fig. El labelled as 
"strict" . 

In order to have again an estimate of the uncertainties of undetermined higher-order contribu- 
tions, we also computed C^ in a thermodynamically consistent manner directly from the pressure, 
given by Eq. (|84|) plus the perturbative zero-temperature result to order 5^^. This has the slight 
deficiency of including higher-order terms in the second term of Eq. (|72|) beyond the accuracy of 
the first one. The "thermodynamically consistent" result is displayed in Fig. [S] by dashed lines, 
and one can see that there is not much difference in the region where C^jC^ is larger than one. 

The results are given for three different couplings. The lines marked "QED" correspond to 
5eflf = 0.303 or uqed ~ 1/137, and the results for g^.^ = 2,3 correspond to as ~ 0.32,0.72 in 
two-flavor QCD. (Recall that g1^ = g'^Nf/2.) While in QED the effect is tiny (the deviations from 
the ideal-gas value have been enlarged by a factor of 20 in Fig. El to make them more visible) , 
in QCD we find that there is an interesting range of T/n where there is a significant excess of 
the specific heat over its ideal-gas value, whereas ordinary perturbation theory 35] would have 
resulted in a low-temperature limit of Cy/C^ = 1 — 2as/iT. 

According to Ref. [13], the critical temperature for the color superconducting phase transition 
may be anywhere between 6 and 60 MeV, so with e.g. a quark chemical potential of ^ = 500 MeV 
the range T/fj, > 0.012 in Fig. El might correspond to normal quark matter. 

Thus, while the effect remains small in QED, it seems conceivable that the anomalous terms in 
the specific heat play a noticeable role in the thermodynamics of proto-neutron stars, in par ticular 
its cooling behavior in its earliest stages before entering color superconductivity J3fll.l4fll.l4ll]. 

If color superconductivity leaves some quark matter components unpaired, the larger values of 
Cj}/C^ at smaller Tf/j, may also be relevant for neutron stars with a quark matter core. 

VII. SUMMARY 

For temperatures much smaller than the chemical potential of quarks (or electrons in the case 
of QED) we have computed the leading contribution to the interaction part of entropy and specific 
heat. For temperatures smaller than the Debye mass oc 5eflF/^) the anomalous (non-Fermi-liquid) 
contributions become dominant. As we have discussed at length, this effect can be viewed either as 
a consequence of a logarithmic singularity of the fermion self-energy at the Fermi surface caused by 
long-range quasi-static magnetic interactions, or more directly as a contribution of the (imaginary 
part of the) transverse gauge boson propagator to the thermodynamic potential when the hard 
fermion degrees of freedom are integrated out first. 

This latter approach proved to be advantageous for a systematic calculation beyond the well- 
known leading-log approximation. We have obtained a hard-dense-loop resummed expression which 
continuously interpolates between the more familiar plasmon effect oc g^ coming from longitudi- 
nal Debye screened gauge bosons and the non-Fermi-liquid effects coming from only dynamically 
screened magnetic interactions. At temperatures much smaller than the Debye mass, we have ob- 
tained a low-temperature expansion starting with the well-known anomalous TlnT behaviour and 
involving also fractional powers of T in subleading terms. The complete HDL-resummed result also 
contains contributions which do not show up at any finite order of the low-temperature expansion, 
being exponentially suppressed by factors of e"'^^' at small T, but which become numerically 
important for intermediate temperatures. 

Finally we have presented a numerical evaluation of our HDL-resummed result together with 
its low-temperature expansion, and we have compared with the exactly solvable large- A'^^ limit of 
QCD and QED. This comparison seems to indicate that our leading-order-in-g( result, which is 
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equally applicable to finite- A'^y QCD, is quite stable in the range of temperature where there is 
an excess of entropy and specific heat over their respective ideal-gas values. In QCD, where the 
coupling as well as the number of gauge bosons is much larger than in QED, the deviation from 
naive perturbation theory is appreciable for T / ^i < 0.05 and thus should be taken into account e.g. 
in studies of thermodynamic properties of quark matter in (proto-)neutron stars. 
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APPENDIX A: CALCULATION OF c 

From the terms which are explicitly shown in Eq. (jSEJ we find the following contribution to the 
coefficient of 6^^ in the pressure, 



7,W 



1 



8107r2 



367r2 + 



(vr^ - 8) ('248 - 967^; - Qvr^ + 48 ln(47r) + 8640^-^ ) . (Al) 



However, some of the terms in the integrand of Eq. H58() that are formally of higher order than 6^^ 
contribute also at the order of 6^^, because the x-integration would be infrared divergent, were it 
not for the cutoff Xmin oc h^. 

Since Xmin depends on g^s only through 6, and since we can drop terms in the integrand of 
Eq. H58|) involving g^^ explicitly, it is sufficient to take the HDL self energy in the following. Then 
the gluon self energy can be written as 

TlT^TlTH-^=9l^^'Hr(^^, (A2) 



q 



„l/3 



with some function Ht- In the following we may neglect the explicit q^ and q^ in Eq. (|53|1 . because 
these two terms do not become singular for small x.'^ After expansion of the integrand with respect 
to h we then obtain integrals of the type 

Now we see clearly that from arbitrary powers of b in the integrand we get contributions to the order 
6^2 in Pt- The case re = 3 corresponds to the term of order b^"^ In b, which we have evaluated already 
in Sect. IIV A[ For n > 3 we can concentrate on the IR region and send to upper integration limit 
in Eq. (|A3|) to infinity. (The cases n < 3 have been evaluated explicitly in Eq. HA1|) .) Furthermore 
we see from Eq. ()A3|) that from the y integration we always get a factor 

^ ^ (A4) 
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However, q and q^ would have to be taken into account when summing up the IR contributions to the coefficient 
of T®, since for this coefficient also less IR singular terms are important [Sg. 



The complete coefficient can thus be written as 
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d" /imn™L\ 

--— arctan ^t;-^ 



(A5) 



fe=0,y=l 



This expression is in fact independent of h (see Eq. ()A,S|) ) . Therefore we may simply set 6 = 1. 
Summing up the (Taylor) series, we find after the substitution x = Anz^ 



(1) ^ f°° f 128 + 371^ z^ -Stt^ (2 + 3z^) 

^^=^^-15/ H ^z 



+z arctan 



vrd - z' 



2z+(z2_l)ln(f±i 



(A6) 



From this expression we see that the complete HDL self energy is required for this coefficient (and 
not only the expansion for small qq, which is sufficient for the fractional powers and the logarith- 
mic terms). The remaining integral over the parameter z can probably not be done analytically. 
Numerically one readily finds 



CT 



-0.00178674305 . , 
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The constant c^ in Eq. ()65() can be determined by summing up IR enhanced contributions in a 
completely analogous manner. The result is 
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~ 0.11902569216... 



(A8) 
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